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Eigenstates of two-dimensional charged electron-hole complexes in mag- 
netic fields are considered. The operator formalism that allows one to par- 
tially separate the center-of-mass motion from internal degrees of freedom 
is presented. The scheme using magnetic translations is developed for cal- 
culating in strong magnetic fields the eigenspectra of negatively charged 
excitons X~ , a bound state of two electrons and one hole. 
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1 Introduction 



Identification of cliarged excitons in magneto-optical spectra of quasi-two- 
dimensional (quasi-2D) systems has induced much interest in the behavior of 
these three-particle electron- hole (e-h) complexes. The negatively, X~, and 
positively, X~^, charged excitons are the bound states of two electrons and 
one hole {2e-h) and two holes and one electron (2/i-e), respectively. In mag- 
netic fields B, in addition to the spin-singlet, higher-lying triplet states of 
X~ and have been observed Theoretically, free charged excitons have 
been studied in strictly 2D systems in the limit of high and low magnetic 
fields and in quasi-2D systems at high magnetic fields . For one-component 
electron systems in magnetic fields, the center-of-mass motion separates from 
internal degrees of freedom. The well-known Kohn theorem which states 
that the electron cyclotron resonance is not shifted or broadened by electron- 
electron interactions, is based on this fact. For e-h systems such a complete 
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separation is not possible in magnetic fields. Nonetheless, any charged in- 
teracting system in a uniform B possesses an exact dynamical symmetry — 
magnetic translations and references therein). It has been shown re- 

cently that due to this symmetry, magneto-optical transitions of charged 
e-h complexes are governed by an exact selection rule, which leads to some 
rather unexpected spectroscopic consequences for charged excitons in B. In 
this work, using an operator formalism, we construct a basis compatible with 
the exact dynamical symmetries — rotations about the 5-axis and magnetic 
translations. Physically, this is equivalent to a partial separation of the center- 
of-mass motion from internal degrees of freedom in B We demonstrate 
that this basis can be used for high-accuracy and rapidly convergent calcula- 
tions of bound X~ states in strong magnetic fields. Our results can also be 
relevant for atomic ions with not too large mass ratios in ultrastrong magnetic 
fields 0. 



2 Basis compatible with magnetic translations 



We consider a strictly 2D system containing two electrons and one hole in a 
perpendicular magnetic field B = (0, 0, B) described by the Hamiltonian 

H = Hq + Hee + Heh , (1) 

HQ=y lk + ^ (2) 

-f^ee = -j T , -f^eh = — ^ -j T , (3) 

e|ri — r2| ^^-^^ 2 ^1^* ~ 

where tcj = —ih'Vj — — A(rj) are kinematic momentum operators. We will 
use the symmetric gauge A = |B x r. The exact eigenstates can be char- 
acterized by the total angular momentum projection M^, an eigenvalue of 
Lz = Z^jlfj X —ift-'^j)z, by the total spin of two electrons Se = (singlet 
states) or S'e = 1 (triplet states), and the spin state of the hole Sh- The latter 
simply factors out and will be disregarded. Performing an orthogonal trans- 
formation of the coordinates {ri, r2, Vh} — > {r, R, r/i}, where r = (ri — r2) / a/2 
is the electron relative and R = (ri + r2)/\/2 center-of-mass coordinates, the 
complete orthonormal basis with a fixed value of can be constructed 
(see also [jlT|) as an expansion in Landau levels (LL's) 



^'il,(r)0£U(R)0S,U(r.). (4) 



Here (pn^ir) = (f)^^{r) are the e- and h- single-particle factored wave functions 
in 5; n is the LL quantum number and m is the oscillator quantum number 
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(see, e.g., For, e.g., zero LL's 



where z = x + is the 2D complex coordinate and Ib = {hc/eBY^'^. The 
factored wave functions are constructed with the help of the oscillator Bose 
ladder operators: For electrons (the charge — e < 0) 

&{r) = ^L={r\{AlnBlr\0) , (6) 
y/nlml 



here the intra-LL operators Bl{rj) = —i^c/2hBe Kj_, where Kj± = Kj^ ± 
iKjy and = ttj — x B (see, e.g., The electron inter-LL operators 

are Al{rj) = —i^c/2hBe7ij^, where '7ij± = njx ± i'^jy The operators com- 
mute as [Ae, Al] = 1, [Be, Bl] = l, and [Ae, Bl] = [Ae, Be] = 0. The analogous 
intra-LL and inter-LL operators for the hole (the charge e > 0) are, respec- 
tively, Bl{rh) = —i^Jc/2hBeKh+ and A\{rh) = —i^c/2hBeTTh-- These can 
be considered as linear functions of spatial coordinates and derivatives and 
have the form 





( z 


dz* 




( 






\21b 



4(r) = 5l(r) = — (;^-2£^^) , (7) 
Bl{v) = Ai{v) 

Single-particle angular momentum projection operators L^e = ^l^e ~ ^iBg 
and Lzh = bIb^ — A^A^, so that m^g = — "n^zh = n — m. The basis (H) includes 
therefore different three-particle 2e-h states such that = ni+r?,2— mi— m2 — 
Uh + ruh is fixed. Permutational symmetry of identical particles requires that 
for electrons in the spin-singlet Sg = (triplet Se = 1) state the relative motion 
angular momentum rii — rrti should be even (odd). The basis (^ proved to be 
effective in strong B for studying impurity-bound states of e-h complexes | |10| ], 



collective excitations — magnetoplasmons and spin- waves |T2|, and effects of 



lateral confinement in quantum dots in B []T3|,|T^ . The equivalent LL expansion 
(using the coordinates {ri,r2,r/i}) has been exploited 0,^] for studying free 
charged excitons in B. However, for translationally invariant systems the basis 
(|^) is not compatible with the magnetic translations. 

Indeed, the Hamiltonian (|T]) commutes with the operator of the magnetic 
translations K = Y.j 0Ji]- Noting that [K^, Ky] = -i^Q, where the total 
charge Q = J2j = ~e for the X~ , one obtains the lowering and raising Bose 
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ladder operators for the whole system 



here k = yc/hB\Q\ K. Therefore, = k+k^ + k^k^ has the discrete oscillator 
eigenvalues 2k + 1, k = 0, 1, . . .. These can be used, together with M^, for 
labelling of exact charged eigenstates of (|1|). Due to the non-commutativity of 
Kx and Ky, there is the macroscopic Landau degeneracy in k. Note now that 

k^ = kj) = J2j k^ + J2i=/=j kj ■ kj is not diagonal in the basis due to 

the cross terms kj ■ kj. 

In order to make the basis (1) compatible with the magnetic translations, 
a canonical transformation diagonalizing k^ should be performed. We deal 
formally with a set of coupled harmonic oscillators. Note then that 

L = Bl{v,) + Bl{v,) - i?,(r,) = V2BliR) - B^^vh) (10) 



and define 

Bl{K)=uBl{R)-vBj,{rh) , B,{K) = uB,{R) ~ vBlir^) , (11) 

where u = ^/2, f = 1. It is this pair of Bose ladder operators in which k^ is 
diagonal: k^ = 2BlBe + 1. Equation ([11]) is in fact a Bogoliubov canonical 
transformation Bl = SBlS"^ generated by the unitary operator (see, e.g.. 



IMm) 

S = exp{e[5e(R)5,,(r^) - Bi(r,05l(R)]} (12) 

with u = ch0 = v^, V = shB = 1. The second pair of linearly independent 
transformed operators -B^(r/i) = SBl{rh)S'' and Bh{rh) = SBh{rh)S^ are 

5j,(r,) = uBlirf,) - vB,{R) , B^irt) = uB^ir^) - vBI{R) . (13) 

The complete orthogonal basis compatible with both axial and translational 
symmetries therefore is 

4(r)"Mt(R)"Mt(r,)"-5t(K)fc5t(r)™5i(r,)'|0). (14) 



In (|14|) the oscillator quantum number is fixed and equals k while = 
—k — m + I + rii + n2 — rih and ni — m is even (odd) for S'g = (S'e = 1). 
The Hamiltonian (0) is block- diagonal in the quantum numbers k,M2,Se- 
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Moreover, due to the Landau degeneracy in k, it is sufficient to consider the k = 
states only. This effectively removes one degree of freedom and corresponds 
to a partial separation of the center-of-mass motion from internal degrees of 
freedom for a charged e-h system in a magnetic field (cf. [^H). 



In (113) the new vacuum |0) = S\0) has been introduced. Disentangling the 
operators in the exponent of S (see, e.g., one obtains 



S = exp {-the BlBl) 

X exp (- ln(che)[S|Se + bIBh + 1]) exp (the B^Bh 



so that 



|0) = S\0) = ^exp {-th&Bl{r,)Bl{R)) |0) 



(15) 



(16) 



For a charged system of N^, electrons and Nh holes (with, e.g., A^'e > Nh), 
a transformation analogous to (|TT])-(|T^ can also be performed. It should 
involve the intra-LL e- and h- center-of-mass operators i?^(Re) and i?/j(R/i) 
with thG = 



Nh/N,; here 



rei/VNe and R/, 



3 X states in lowest Landau levels 



We now demonstrate how the developed formalism works. We will consider 
the limit of high magnetic fields P,^,^ 



TT e 
27h 



(17) 



when mixing between different LL's can be neglected. Eq is the characteris- 
tic energy of the Coulomb interactions in strong B, TiuJcc{h) = heB/me(h)C. 
Charged magnetoexcitons can then be labeled by the total electron LL num- 
ber He = ni + n2 and by the hole LL number nh. Indeed, when (17) is fulfilled, 
the states having different quantum numbers UeUh and n'^n'f^ are only weakly 
~ £"0/1 (n'g — ne)huce + {n'h — nhjhu ch\ mixed by the Coulomb interactions ||TT|| . 



We focus on the states in zero LL's [ni = n2 = = Q m. ([T^) ]. The operators 
p!3D have a simple representation in the new coordinates Pi = \/2 R — 
V2rh - R: 5t(R) = Bl{p,) and Blirj,) = Bl{p,). The complete 

and arbitrary Mz = 



m 

and p2 

infinite orthonormal basis in zero LL's with fixed k 
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/ — m takes the form 

^J^Bl{rrBi{p,y\0)^\ml) (18) 

with odd m = 2p+l (even m = 2p), p = 0, 1, ... in the electron triplet 3^ = 1 
(singlet Se = 0) states. The Coulomb interactions in the new variables are 

Hee = /- , -f^eh = T p— ^ T ■ (19) 

The matrix elements of the e-e interaction are diagonal in the basis ([T8|) : 

where Vo,m is the interaction of the electron with a fixed negative charge — e in 
zero LL (e.g., |]Tn[)- Due to the permutational symmetry, the two terms in H^h 
give the same contributions; calculations, however, are not so straightforward 
as (PO]). This is connected with the fact that the coordinate transformation 
{r, R, r/i} —>■ {r,p^,p2} is not orthogonal. As a result, the coordinate repre- 
sentation of the new vacuum is not factored in and pg- 

(rPiP2|0) = ,,,,,, exp -2 , (21) 



v^(27r£|)3/2 



here Zj = pj^ + ipjy, j = 1, 2. Therefore, when acting on the vacuum |0), the 
operator Bl{p2) = B\^{Yh) gives a combination {Z2 + -^Zi) / \/2^B = z^IUb- 
To eliminate the coordinate p^, we perform the shift p^ ^ p = p^ + ^p2 = 



■^R and obtain 



,2 



2£|, 




Zi-miJ2-ni2 • 



Integrating out the variable p, we reduce the problem to an effective two- 
particle e-h problem in zero LL's (cf. |1^). The peculiarity of the situation 
is that the effective particles are characterized by different magnetic lengths. 
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The matrix elements ( p2D can be presented in the form (mi = m, m2 = m + s, 

h = I, I2 = I + s) 



[m 



+ sl + s\H^^\ml) = (-2V2)2-'-f ^ (C^C^::)" ut='\s) , (23) 



fe=0 



where are binomial coefficients and the matrix elements of the Coulomb 
interparticle interactions in zero LL's have been introduced: 



rfVi / d\2 00^2 (ri) 0o£* (1^2/ Va) - — r0S,fcJ (rs/ V") (1 

J 1 1"! — T2\ 

(a) 



5fci-mi,fe2-m2f^min(fcl,fc2),min(mi,m2)(l"^l —"^2!) (24) 

The matrix elements (p^ ) can be found analytically for arbitrary a: 



(1 + a)"+2 2™+"+^ fc^o «=o 



^(IT^ [2(fc + / + .)- 1]!! [2(m -k)- 1]!! [2(n - /) - 1]!! 



Equations (pO]), (|23|), and (pSj) determine the secular equation of the infinite 
order that should be solved to obtain the three-particle 2e-h states in zero 
LL's. A truncation of the basis should naturally be performed. An important 
property of the developed basis (0) [and ([T8|) ] is that such a truncation does 
not break the translational invariance. On the contrary, a truncation of the 
basis (Q), as performed in (see also p|), leads to spurious mixing of 

different /c-states and violates the exact magneto-optical selection rule — 
the conservation of the oscillator quantum number k. 

The developed approach also provides an effective computational tool: First, 
we have been able to remove one degree of freedom in the three-particle prob- 
lem, so that configurational space is substantially reduced (cf. p). As a result, 
with finite-size calculations it is even possible to reproduce with a reasonable 
accuracy the three-particle continuum — a neutral magnetoexciton plus a 
scattered electron |Q. Second, for bound X~ states lying outside the continua 
we have extremely rapid convergence within each LL. This is associated with 
the exponential decay of the off-diagonal matrix elements (^3]). Consider, e.g., 
the k = triplet Xj~^^q„^^o state in zero LL's with = —1. The asymptotic 
behavior of the relevant off-diagonal Coulomb e-h matrix elements is 



2v^, (a=2)/„ N / 32 /I 



(2s + 12s|i/eh|10) = -^f/^r^^(2s)^-W— (^-) i?o,s»l. (26) 
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Table 1 

Singlet X~^^^^ and triplet X^^^^ charged magnetoexcitons in Landau levels ngn/i 





M, 


Interaction energy {Eq) 


Binding energy (Eq) 




-It 


-1.04345 


0.04345 




-3t 


-0.78056 


0.20690 




-4t 


-0.75776 


0.18410 




It 


-1.08596 


0.08596 



t the only bound states in the given LL's neUh. 

t the ground states among many other bound states with the same neiih. 

Also, even the 1x1 matrix Hamiltonian in the basis (|18D (1 0|iJcc + -f^ch|l 0) = 
— 1.0073^^0 ensures the X~ binding: it gives a positive binding energy 0.0073i?o; 
this is relative to the ground state energy —Eq of the neutral X„^=onh=o mag- 
netoexciton in zero LL's. As a result, the X^q binding energy can be calculated 
with virtually unlimited accuracy and equals 0.043452ii^oi this value is com- 
patible with 0J^. Not accounting for the Landau degeneracy in k, the X^q 
state with = —1 is the only low-lying bound X~ state in zero LL's: there 
are no other bound triplet or sing let states ||,|J|]. 

Similar considerations apply to the X~ states in higher LL's. Some of the 
results for the X~ ground states are presented in Table 1. There is only one 
bound X~ state in the first electron LL (the basis (1T3) includes the states 



with rii = 1, = 0, = and rii = 0, ?t,2 = 1, Uh = 0). This state 
is the triplet X^^o with Mz = 1, whose binding energy is almost twice that 
of the Xjoo state in zero LL's 0. This resembles a stronger binding of the 
triplet D" state (two electrons bound by a donor ion) in the first electron 
LL [|n[] and has the same physical origin. The X^^o binding energy is counted 
from the lowest possible unbound state in the same LL's, which is the neutral 
magnetoexciton X„^=onh=o with the second electron in the scattering state 
in the rig = 1 LL. As calculations show, there are many bound X^ states 
in the next hole LL [rii = n2 = 0, n/i = 1 in (|T^)] — both triplets X^gi and 



singlets X^^ (see Fig. 1). These are lying below the ground state of the neutral 
magnetoexciton X„^=onh=i; the latter has the energy — 0.57366i?o- Due to this 
small binding energy of the neutral Xn^=on,,=i magnetoexciton (comparatively 
to the X„^=onh=o magnetoexciton), the triplet X^q^ and singlet X^^^ ground 
states have rather large binding energies (Table 1). In all LL's, there are also 
higher-lying bound three-particle 2e-h states originating from the internal 
bound motion of 2D electrons in strong magnetic fields. These states appear in 
the spectrum at relatively large positive values of the total M^, that correspond 
to the hole being at large distances from the electrons (cf. with the similar 
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Fig. 1. Bound and scattering three-particle 2e-h singlet Se = states in the 
(rae"h)=(01) LL's. Note that only the k = states are shown. The energy is given 
relative to the energy of free LL's ^hujce + |^Wch, in units of £"0 = \/7r/2 e^/e/e- 
Open circles correspond to low-lying X~q^ states. Filled circles correspond to ex- 
cited states originating from bound internal motion of two electrons in 2D in high 
magnetic fields. The shaded area shows the continuum corresponding to the neutral 
magnetoexciton Xqi plus a scattered electron in the rig = LL; the lower continuum 
edge lies at an energy — 0.5737£'o. 



states in the D problem ITT 



4 Summary 



In conclusion, we have developed a formalism that allows one to preserve the 
exact symmetry — magnetic translations — when performing the Landau 
level expansion for charged electron-hole complexes in magnetic fields. This is 
achieved by using the Bogoliubov canonical transformation mixing the center- 
of-mass motions of the electron and hole subsystems. The effectiveness of 
the scheme has been demonstrated for high-accuracy and rapidly convergent 
calculations of two-dimensional charged excitons X~ in magnetic fields. This 
can be useful for studying the eigenspectra of charged excitons in quasi-two- 
dimensional quantum wells at strong and intermediate magnetic fields. 
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